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We perform a comprehensive statistical analysis of the standard model (SM) with three and four 
generations using the latest Higgs search results from LHC and Tevatron, the electroweak precision 
observables (EWPOs) measured at LEP and SLD and the latest determinations of Mw, m t and a s . 
For the three-generation case we analyse the tensions in the electroweak fit by removing individual 
observables from the fit and comparing their predicted values with the measured ones. In particular, 
we discuss the impact of the Higgs search results on the deviations of the EWPOs from their best- 

+3.0 

fit values. Our indirect prediction of the top mass is m< = 175.7_ 2 2 GeV at 68.3% CL, in good 
agreement with the direct measurement. We also plot the preferred area in the Mw-m t plane. The 
best-fit Higgs mass is 126.0 GeV. For the case of the SM with a perturbative sequential fourth 
fermion generation (SM4) we discuss the deviations of the Higgs signal strengths from their best-fit 
values. The H — > 77 signal strength now disagrees with its best-fit SM4 value at more than 4a. We 
perform a likelihood-ratio test to compare the SM and SM4 and show that the SM4 is excluded at 
5.3cr. Without the Tevatron data on H — > bb the significance drops to 4.8a". 



INTRODUCTION 

Electroweak precision observables (EWPOs) played an 
important role in the prediction of the mass of the top 
quark prior to its discovery [TJH]- Later, with improving 
Tevatron data on the mass m t of the top quark, EW- 
POs were used to constrain the mass ijih of the Higgs 
boson, albeit with little precision since EWPOs only de- 
pend logarithmically on this quantity Recently AT- 
LAS and CMS have discovered a convincing candidate 
for the Higgs boson with a mass around 126 GeV [H[S]. 
With the information on ran available, EWPOs enter a 
new era, as they directly test the standard model (SM) 
without involving otherwise undetermined fundamental 
parameters. In this paper we present a combined fit of 
the EWPOs and Higgs signal strengths in the decays to 
77, WW, ZZ, bb, and rr studied at the LHC and the 
pp — » H — > bb signal strength determined at the Tevatron 

!• 

The SM is minimal in several respects, e.g. the fermions 
belong to the smallest possible representations of the 
gauge groups and electroweak symmetry breaking is 
achieved with a single Higgs doublet. However, the 
fermion field content is non-minimal and organised in at 
least three families. There has been a tremendous inter- 
est in the phenomenology of a fourth fermion generation, 
with more than 500 papers in the last decade. The SM 
with a sequential fourth generation, SM4, had survived 
global analyses of EWPOs and flavour observables [7HT5] , 



but was put under serious pressure from the first LHC 
data on Higgs searches [THll8| . In this paper we show 
that the perturbative SM4 is the first popular model of 
new physics which is ruled out by the LHC at the 5a 
level. This strong statement is possible because of the 
non-decoupling property of the SM4, as loops with 4th- 
generation fermions do not vanish with increasing masses. 
For the same reason it is difficult to compute the statis- 
tical significance at which the SM4 is ruled out: the non- 
decoupling property implies that the SM4 and SM are 
non-nested, i.e. the SM is not obtained from the SM4 by 
fixing the additional parameters. This complicates the 
statistical procedure which compares the performance of 
the two models at describing the data. We have solved 
this problem with the help of a new method for toy Monte 
Carlo simulations, as implemented in the mj/Fitter pack- 
age [in] ■ A first application of mj/Fitter to the SM4 with 
the data available before the announcement of the Higgs 
discovery right before the ICHEP 2012 conference has 
been presented in [TB] . 



METHOD AND INPUTS 

We combine electroweak precision data with Higgs sig- 
nal strengths provided by ATLAS, CMS and Tevatron. 
(To study the impact of the excess in H — > bb events 
reported by CDF we also show results with the CDF in- 
put excluded.) Our fit parameters in the SM are the 
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TABLE I. Experimental inputs for Higgs signal strengths. 
Except for — > 77 CMS only provides signal strengths at 
125.5 GeV. ATLAS has not published a 2012 update on 
H — > tt and H — > bb, so we take the 2011 data. 

Z mass Mz, the top quark mass m t , the strong cou- 

(5) 

pling a s , the hadronic contribution Aa^ to the fine 
structure constant at the scale Mz in the five-flavour- 
scheme and the Higgs mass run- In the SM4 the addi- 
tional parameters are the fourth generation quark masses 
mf and my and the fourth generation lepton masses 
mg i and m„ 4 . Mixing between the quarks of the fourth 
and the first three generations is neglected in this anal- 
ysis, since our previous study |18] showed that such 
mixing is disfavoured by the combination of Higgs sig- 
nal strengths and EWPOs. Direct searches for fourth- 
generation quarks at the Tevatron and the LHC put lower 
bounds on the masses of the fourth-generation quarks. 
For example, the current highest limit on rrib> is given 
in [2U] as my > 611 GeV. However, these limits rely 
on specific assumptions about the mass splitting and the 
decay patterns of the heavy quarks. We therefore use 
a conservative limit of m t ' , my > 400 GeV in our fits. 
For the lepton masses we require 771^4 > 100 GeV and 
to„4 > Mz/2. As upper limit for all fourth generation 
fermion masses we choose 800 GeV. 

Our inputs for the Higgs signal strengths are sum- 
marised in Table [T] For a given Higgs production and de- 
cay mode X — >• H —> Y, the signal strength jl(X — > H — » 
Y) is defined as the observed production cross section 
times branching ratio divided by the SM prediction. The 
asymmetric errors are accounted for by using an asym- 
metric gaussian likelihood function. For the SM fit and 
the Higgs mass scans we treat mu as a free parameter 
and interpolate the data from signal strength plots ver- 
sus Higgs mass, as provided by the ATLAS, CMS, CDF 
and DO collaborations. When comparing the SM and the 
SM4 we keep the Higgs mass fixed in our fit and use the 
combined signal strengths given in Table [T] as inputs. 

On the theory side, the SM4 signal strengths are com- 
puted by appropriately scaling the SM branching frac- 
tions and production cross sections separately for each 
production mechanism. (Further details can be found in 
|16|.) In this sense, our treatment of the Higgs signal 
strengths is a special case of effective coupling analyses 



such as [251427J . Unfortunately, these analyses are in- 
sufficient for constraining, let alone ruling out a concrete 
model like the SM4. An effective coupling analysis which 
"contains" the SM4 would have to treat the Higgs cou- 
plings to 77, WW, ZZ, gg, bb, tt and v±v± as indepen- 
dent parameters and provide full information about the 
X 2 function on this seven-dimensional parameter space. 
Even then one could not compute the p- value of the like- 
lihood ratio test comparing the SM and the SM4, since 
this requires a numerical simulation with toy measure- 
ments. Hence the results of [23 -27 cannot be applied to 
the SM4. 

Although our experimental inputs are combined results 
from the 7 and 8 TeV LHC runs, we compute the signal 
strengths using 7 TeV SM cross sections only. This prac- 
tice is justified because the signal strengths only depend 
on the ratios of Higgs production cross sections for dif- 
ferent production mechanisms and not on their absolute 
size. The ratios are constant to a good approximation 
when going from 7 to 8 TeV [28 . Note however, that we 
treat the H — > bb signal strengths from the Tevatron and 
LHC detectors as two different observables because the 
ratio of W and Z associated production cross sections is 
different at the Tevatron and at the LHC. 

Heavy 4th generation fermions imply large Yukawa 
couplings which eventually make the theory non- 
perturbative. The 1978 paper [2S] estimated a break- 
down of perturbation theory at my > 500 — 600 GeV 
from considerations of tree-level partial wave unitar- 
ity |30j . However, this bound merely implies that for 
my « 500 GeV loop corrections become important. In 
our fits we compute the Higgs width and branching ratios 
in the SM4 and SM with HDECAY v. 4.45 [51], which 
implements the higher-order corrections of [52TT55] (see 
also [56]). 

The global fits with a variable Higgs mass were done 
with the CKMfitter software [37]. The EWPOs in the 
SM4 were calculated with the method described in |38) . 
using FeynArts, FormCalc and LoopTools [39 -WT\ to 
compute the SM4 corrections to the EWPOs. The EW- 
POs in the SM were calculated with the ZFitter soft- 
ware [42TI44] . The SM Higgs production cross sections 
were taken from [45] (LHC) and [46] [47] (Tevatron). For 
the numerical computation of the p-values we use the 
mj/Fitter package [TH] which in turn uses the Dvegas code 
[4"B"H5n] for numerical Monte Carlo integration. 

SM FIT RESULTS 

The EWPO inputs as well as the SM fit results can be 
found in Table [IlJ The experimental error on the Higgs 
mass is determined in a fit to all available signal strengths 
(see the second column of Table |T| . In addition to the 
best-fit value and 68.3% CL interval of each observable 
we also show the prediction for each observable, which 
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TABLE II. Experimental inputs and fit results for the electroweak precision observables in the SM. The observables were 
calculated with ZFitter |42H44| . The inputs are listed in the second column. The first error is statistical while the second (if 
present) is systematic. We also use the correlations from [52]. The input for a s is the value determined from the r lifetime. 
In the fourth column, we show the results of a global fit using all available inputs. Here, the errors are 68.3% CL intervals. 
The fifth column contains the prediction for each observable, obtained by removing the direct input for that observable and 
re-running the fit. The corresponding difference of the minimal \ 2 values is shown in column six. The quantities in the last 
five rows were used as fit parameters. 



is obtained by removing the direct input for that observ- 
able and minimising the \ 2 function again. The resulting 
change Ax 2 in the minimum \ 2 value is shown in the 
last column of Table [ilj A large A\ 2 value indicates that 
the observable is in strong disagreement with the other 
observables. The most prominent "outlier" is therefore 
Ap B , followed by cr° ad and A\. The leptonic left-right 
asymmetry Ai is the main cause for the disagreement 
between the predicted and measured value of mn- re- 
moving the Ai input leads to a predicted Higgs mass of 
124 GeV. Note however, that the updated inputs for to* 
and M\y move the predicted Higgs mass up to 108 GeV. 

The deviations of the EWPOs from their best-fit pre- 
dictions are shown in Figure [T] This figure shows the fit 
with Higgs signal strength and EWPO inputs as well as 
the fit with EWPO inputs only. Note that, due to the log- 
arithmic dependence of the EWPOs on the Higgs mass, 
the inclusion of the Higgs signal strengths is essentially 
equivalent to fixing the Higgs mass at 126 GeV. We see 
that the new Higgs data has a relatively small impact on 
the deviations of most EWPOs. The main difference is 
an increase in the deviation of M\y to 0.8cr. The 68.3%, 
95.5% and 99.7% CL regions in the mt-Mw plane (using 
Higgs signal strengths and EWPOs) are shown in Figure 




SM4 FIT RESULTS 

The impact of a fourth fermion generation on the Higgs 
signal strengths has been discussed extensively in the lit- 
erature. The Higgs production cross section via gluon 
fusion is enhanced by a factor of 9 due to the additional 
heavy quarks in the loop 0[S!|]- In H — > 77 searches, this 
factor is overcompensated by a reduction of the branch- 
ing ratios, which is due to an accidental cancellation be- 
tween gauge boson and fermion loops at next-to-leading- 
order [33] . Finally, all signal strengths can be suppressed 
by a common factor if the invisible H — > 1/41/4 decay is 
kinematically allowed [60Tf67]. 

The deviations of the Higgs signal strengths in the SM 
and the SM4 are shown in Figure |3j For comparison 
we also show the deviations for the SM4 fit to the data 
available before the announcement of the Higgs discovery, 
which was used in [18]. For the SM4 fit to current data 
we also show the shifts in the minimum x 2 value obtained 
by removing individual signal strengths from the fit. We 
see that the deviation of the H —¥ 77 signal strength 
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FIG. 1. Deviations of the EWPOs in the Standard Model. 
The observables were calculated with ZFitter |42H44) . For an 
observable O with experimental value O cxp , experimental er- 
ror AO CX p and best-fit prediction Om we define the deviation 
as (O e xp - Ofit)/AO Hp . 



has increased dramatically with the new data and now 
exceeds four standard deviations. Furthermore, the SM4 
cannot explain an excess in H — » bb searches because 
the Higgs production mechanisms for these searches are 
HW and HZ associated production, which are not en- 
hanced by a factor of 9 like the gluon fusion production 
mode. Thus the fit improves significantly if the Tevatron 
measurement of the H — > bb signal strength is removed. 

Figure [4] shows the minimum x 2 values in the SM and 
the SM4 as functions of the Higgs mass. The absolute 
minimum in the SM4 is at m# = 124.5 GeV and the 
minimum % 2 value is larger than the one in the SM by 
20 units. 

To compute the statistical significance at which the 
SM4 is ruled out one has to perform a likelihood-ratio 
test. This task is complicated by the fact that the SM and 
the SM4 are not nested, i.e. the extra parameters in the 
SM4 cannot be fixed in such a way that all observables 
assume their SM values. As explained in [19], analytical 
formulae for p-values are not valid in this case and one 
has to rely on numerical simulations. In our analysis 
we used the improved simulation methods implemented 
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FIG. 2. The 68.3%, 95.5% and 99.7% CL regions in the m t - 
Mw plane using Higgs signal strengths and EWPOs. Also 
shown are the experimental values of m t and Mw and their 
errors. The inner error bars are the statistical errors. 
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FIG. 3. Deviations (defined as in Fig. [T| of the Higgs signal 
strengths for the SM (blue) and for the SM4 (red) at a fixed 
Higgs mass of 126 GeV. For comparison the results of the fit 
to pre-ICHEP2012 data from [18] are also shown in green. In 
the right column we show, for the SM4 fit to current data, 
the change in the minimum x 2 value when the corresponding 
signal strength is removed from the fit. 



in the mj/Fitter package. For performance reasons, we 
fixed the SM parameters Mz, m t , a s , Aa^ ad and ran 
to their best-fit values in these simulations. This is a 
valid approximation since the SM4 fit is now dominated 
by the Higgs signal strengths and their dependence on 
the SM parameters is negligible. 



summarises 



Table [Hi 

the results of the likelihood-ratio tests. If all inputs are 
used, the SM4 is excluded at 5.3 standard deviations. 
If the Tevatron input for the H — > bb signal strength 
is removed the number of standard deviations drops to 
4.8. Note that these significances hold for an SM4 with a 
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FIG. 4. Higgs mass scan for the SM (blue line) and the SM4 
(red line) based on the input set in the second column of Table 
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lytical formulae for p-values do not hold in the case of 
non- nested models. Using a conservative lower limit of 
400 GeV for the fourth-generation quark masses and fix- 
ing the SM parameters to their best-fit values we find 
that the SM4 with a minimal Higgs sector is ruled out at 
5.3cr. If the results of H — » bb searches at Tevatron are 
excluded from the analyses the SM4 is still ruled out at 
4.8cr. 
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inputs 



p-value standard dev. 



all inputs 1.1 • 10~ 7 

H -> bb from LHC only 1.9 • 1(T 6 



5.3 
4.8 



TABLE III. The p- values and number of standard deviations 
of likelihood-ratio tests comparing the SM and the SM4. The 
SM parameters were fixed to their best-fit values in these 
simulations. 



minimal Higgs sector and may be weakened if the Higgs 
sector of the SM4 is extended 



CONCLUSIONS 

We performed a combined fit of the parameters of the 
standard model with three and four generations, com- 
bining Higgs search results and electroweak precision 
data. In the SM electroweak fit the prediction for the 
Higgs mass from EWPOs has moved closer to the value 
favoured by direct Higgs searches due to new inputs for 
m t and Mw When the Higgs signal strength inputs are 
combined with the EWPOs the discrepancy between the 
measurement of Mw and its best fit value increases but 
stays below la. All other deviations are essentially un- 
affected by the new input. 

In the SM4 the measured H — »• 77 signal strength dis- 
agrees with the best-fit prediction by more than four 
standard deviations. Another source of tension is the 
excess in H —> bb searches at the Tevatron in combina- 
tion with the deficit in H — » tt events. The dominant 
Higgs production mechanism for H — > tt searches (gluon 
fusion) is enhanced by a factor of 9 in the SM4 while 
the relevant production mechanism for H — > bb searches 
(HW and HZ associated production) is slightly reduced. 
The statistical significance at which the SM4 is excluded 
must be computed by numerical simulation methods like 
those implemented in the myFitter software, since ana- 
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